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ABSTRACT 


A  five- level  baroclinic  model  for  numerical  prediction  of  height 
fields  at  1000*  850*  700,  500,  and  300  mb  is  developed  and  programmed. 
Several  forms  of  the  model  are  tested  to  determine  the  contribution 
of  the  various  terms  of  the  vorticity  equation.  Forecasts  from  three 
different  methods  are  compared.  The  model  shows  utility  in  prognosis. 
Although  there  are  individual  differences  among  the  methods,  no  one  of 
the  three  is  significantly  better  over  the  entire  hemisphere  on  the 
data  tested  here. 
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1.  Introduction, 


Until  very  recently*  the  baro tropic  models  have  been  the  principal 
source  of  numerical  weather  prognoses,  Their  shortcomings  in  areas  of 
baroclinic  development*  plus  the  advent  of  sophisticated  digital  computer 
systems*  have  spurred  interest  in  baroclinic  models.  It  is  the  purpose 
of  this  research  to  develop  and  program  a  five-level  baroclinic  model  and 
to  investigate  the  contribution  of  the  various  terms  of  the  vorticity 
equation  to  the  prognosis  of  the  pressure  field* 
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2.  Background, 

In  the  model  developed  here  *  prognosis  is  accomplished  by  means  of 
the  vorticity  equation >  which ,  with  pressure  as  the  vertical  coordinate  * 
can  be  written 


H  =-j  w) 


(M, 


a* 


At 

7>r 


I 


Ml 


sp 


5 /  37  3F 


3fr> 

a* 


JLA-) 


(1) 


Basic  to  the  derivation  of  this  equation  is  the  Helmholtz  theorem  which 
expresses  the  wind  as  the  sum  of  a  divergent  and  a  non- divergent  compon¬ 
ent; 

V4  ^  Vf  +  V^.  (2> 

The  potential  function  ~X^  for  the  divergent  part  of  the  wind  is  obtained 
through  solution  of  the  continuity  equation 

VaX  +  I (3) 

The  stream  function  for  the  non- divergent  portion  is  approximated  by 
solution  of  the  balance  equation 


as  described  by  G,  'Arnason  [i]a  where  (jj>  is  the  geopotential.  Use  of 
the  stream  function  in  the  vorticity  equation  enables  one  to  avoid  the 
geos  trophic  approximation  and  the  concomitant  spurious  anticyclogenesis. 
The  vertical  component  of  relative  vorticity*  X  *  is  the  Laplacian 
of  the  stream  function 

l-V'Y  <5 
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and  absolute  vortlcity,  '[  ,  Is  the  sum  of  this  and  the  vertical  component 
of  the  earth's  vortlcity,  f. 

#  *  r  +  -f  (6) 

Vertical  velocity,  60,  is  computed  by  the  diagnostic  equation 


developed  by  Clarke  and  Lawniczak  [ 3 J  .  The  forcing  function  is  calculated 
with  the  geostrophic  approximation.  This  results  in  some  inconsistencies 
in  the  model  but  to  avoid  reprogramming  the  entire  equation,  we  have  gone 
ahead  with  comparison- type  experiments.  Vertical  velocity  is  assumed  zero 
at  100  mb,  but  other  boundary  conditions  are  among  the  options  available 
to  the  user  of  this  computer  program  as  follows; 

a)  the  lower  boundary  can  be  level  or  include  terrain; 

b)  vertical  velocity  can  be  set  equal  to  zero  at  the  lower  boundary, 
or  computed  with  or  without  terrain  effects  and  with  or  without  the 
effect  of  friction; 

c)  the  stability  parameter  CT  can  be  a  point  variable,  a  function  of 
pressure  alone,  or  a  constant;  and 

d)  the  maximum  acceptable  residual  in  the  3-dimensional  relaxation 
can  be  selected  by  the  user. 

In  the  vortlcity  equation,  the  contributions  from  the  twisting  term, 

/  Ioj  c>^  I'Y  \ 

[  T"*”  c>x  a  p  — Ty  SVd p  )  anc*  t^ie  vertical  advection  of  vortlcity  term, 

UJ  are  usually  considered  small,  and  some  investigations  have 

e)  p  1 

shown  that  they  tend  to  compensate  each  other  [2],  Also,  integral  theo¬ 
rems  show  [4]  that  to  avoid  the  mean  generation  of  vortlcity  in  a  closed 
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volume,  the  two  terms  should  be  included  or  neglected  together.  Both 
have  been  omitted  from  the  initial  investigation  reported  here, 

A  further  simplification  involves  omission  of  the  advection  of  abso¬ 
lute  vorticity  by  the  divergent  portion  of  the  wind,  * 

For  consistency,  it  is  necessary  then  also  to  substitute  the  mean  value 
of  the  Coriolis  parameter,  f  ,  for  the  absolute  vorticity  in  the  diver- 
gence  term,  f  . 

In  the  equations  that  follow,  the  stream  function  has  been  redefined 
so  as  to  be  in  units  of  height, 

*Vr  ^y-E  .  <8) 

The  resulting  forms  of  the  vorticity  equation,  expressed  in  finite- 
difference  notation,  are 


rm1  Q 

"dry 


v 


i  AE-  9  frr? 
At  “  {Hd1 


jj(L>?)+r 


A  OJ 

Ap, 


(9) 


9  tt  / r  yA  q /mV AT-  AT?  ^  & u> 

gp  ^  y  At  |  Jj  fU  p  Lf  \  a*  Ay  AY /  'l  a  p  , 


(10) 


and  the  complete  equation 


rrnx  9  rr?^  aE-  .  -  3/W  Tf/C  t,-,  ^  _  3/W-/  A~XA^  AXA>?\  a,  A  A) 

'  d*  T  *  At  f  4  c!1  ^  ** '  f  t  'Av  AYAV/  :  Ap 

0  9  /V*7  A^  \ 

f  d1  AP  f  4  ^AY  AV  Ap  A*  AYAp'j 

where  d  is  the  grid  distance  and  /W\  is  the  ratio  of  this  distance  to  the 
corresponding  geographical  distance. 
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Horizontal  spatial  operators  are  approximated  by  the  usual  centered 


differences  *  To  approximate  the  vertical  derivatives  in  the  vorticity 
equation  a  Lagrangian  parabola  was  fitted  through  three  points.  For  a 
centered  difference  the  equation  is 


j.  (P.-PO 


(12) 


Ap  p-.p,  CP-P.Hft-Ri) 


All  computations  were  done  on  the  Control  Data  Corporation  model 
1604  computer,  using  the  U.  $.  Navy  Fleet  Numerical  Weather  Facility 
octagonal  grid.  This  grid,  centered  at  the  pole  and  inscribed  in  the  ION 
latitude  circle,  consists  of  1977  grid  points  with  a  grid  size  of  381  km 
at  SON. 

The  week  of  22  November  1962  was  a  period  of  considerable  cyclogene¬ 
sis  in  dense-data  areas,  and  therefore  excellent  for  testing  a  baroclinic 
model.  Initial  data  used  were  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  height  analyses  at  1000,  850,  700,  500,  and  300  mb  for  12002  on 
22  November  1962*  The  height  analyses  for  Q000Z  and  1200Z  23  November 
were  used  in  the  evaluation  of  the  predictions. 
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3, 


Procedures . 


In  the  diagnostic  63 -equation  for  the  initial  tests  discussed  here, 
the  boundary  conditions  used  are: 

a)  vertical  velocity  equal  to  zero  at  100  mb; 

b)  vertical  velocity  equal  to  zero  at  the  lower  boundary; 

c)  no  terrain  effects  at  the  lower  boundary;  and 

d)  the  stability  parameter,  or  ,  is  a  function  of  pressure.  The  pro¬ 
gram  is  used  in  "forecast  mode";  that  is >  temperatures  are  computed  hydro¬ 
statically,  rather  than  taken  from  objective  analysis.  The  maximum 
residual  in  the  3-dlmensional  relaxation  is  1  x  10"^  rab  per  second, 

The  sequence  of  calculations  programmed  can  best  be  explained  with 
reference  to  figure  1.  Input  data  in  every  case  are  the  U,  S*  Navy  Fleet 
Numerical  Weather  Facility  height  analyses  at  the  five  standard  levels, 
modified  to  ensure  vertical  consistency. 

In  the  first  method  with  equation  (9),  stream  functions  are  used 
throughout.  The  inversion  to  height  fields  [l] ,  effectively  the  reverse 
of  the  balance  equation,  is  thus  bypassed  and  stream  functions  are  used 
for  input  data  for  both  diagnostic  and  prognostic  equation.  This  is 
hereafter  referred  to  as  method  I,  The  resulting  24-hour  prognoses  are 
displayed  on  charts  6  through  10  and  discussed  in  section  4  below. 

Since  the  geos trophic  approximation  is  used  in  the  forcing  function 
of  the  63 -equation,  and  temperatures  are  computed  hydrostatically,  there 
might  be  some  reason  to  suspect  the  vertical  velocities  computed  from 
height  fields  could  be  more  accurate  than  those  from  stream  functions-  A 
completely  geos trophic  model  was  programmed,  using  pressure-height  data 
throughout.  This  method  showed  the  expected  spurious  an ticyclo genesis 
even  at  12  hours  and  was  carried  no  further. 
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The  complete  process  diagrammed  in  figure  1  has  been  used  in  the  re¬ 
mainder  of  this  investigation.  In  the  first  time  step  the  initial  data 
(five  vertically-consistent  height  fields)  are  used  as  Input  to  the  bal¬ 
ance  equation  and  to  the  CO -equation.  From  the  former  are  obtained  stream 
functions  for  the  five  levels;  the  latter  yields  vertical  velocities  at 
the  same  levels.  With  these  stream  functions  and  vertical  velocities  the 

vorticity  equation  (equation  (9)  for  method  II,  equation  (10)  for  method 
£  E 

III)  gives  -  ,  the  local  change  in  stream  function,  at  the  five  levelSc 

At 

When  these  are  added  to  the  initial  values,  the  first  forecast  time  step 
is  complete.  As  shown  in  figure  1,  these  new  E  values  are  now  input 
fields  for  both  the  vorticity  equation  and  the  balance  equation;  the 
latter  being  used  for  inversion  to  pressure  height  as  required  by  the 
CJ -equation. 

Subsequent  time  steps  follow  the  same  routine,  except  in  the  applica¬ 
tion  of  A £  ,  which  is  done  according  to  the  centered-dif ference  equation 


At 


(13) 


In  order  to  minimize  the  errors  inherent  in  the  initial  forward 
time  step,  the  first  time  increment  is  only  15  minutes.  Centered  time 
steps  are  then  used  to  complete  successive  15-minute  prognoses  to  one 
hour,  and  finally  1-hour  prognoses  to  12  hours.  The  fields  are  then 
smoothed  lightly  to  remove  any  small  perturbations  resulting  from  trunca¬ 
tion  or  round-off  errors,  and  the  process  is  repeated  to  24  hours. 
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4, 


Results . 


In  review;  method  I  uses  stream  functions  throughout  *  and  with 
equation  (9)  as  the  prognostic  equation  excludes  the  advection  of  abso¬ 
lute  vorticity  by  the  divergent  portion  of  the  wind  and  substitutes  f 
for  7l  in  the  divergence  term;  method  II  inverts  the  stream* function 
fields  to  pressure  heights  at  every  time  step,  but  retains  the  other 
simplifications  of  method  I;  method  III  also  inverts  to  pressure  height 
at  each  time  step,  but  uses  the  more  complete  equation  (10)  as  the  prog* 
nos tic  equation. 

Charts  1-5  are  the  pressure  height  analyses  at  the  beginning  of  the 
period.  The  next  15  charts  are  the  results  of  the  three  24*hour  fore* 
casts;  method  I  in  charts  6-10;  method  II,  charts  11-15;  and  method  III, 
charts  16-20,  Charts  21-25  are  the  observed  height  analyses  at  the  end 
of  the  forecast  period.  Chart  26  is  the  U,  S.  Navy  Fleet  Numerical 
Ueather  Facility  baro tropic  forecast  at  500  mb  for  the  same  period. 

Table  1  shows  the  contour  origin  and  interval  for  all  the  charts. 

In  general,  corresponding  charts  from  the  three  methods  are  quite 
similar.  They  all  retain  a  typical  meteorological  appearance  and  the 
various  systems  show  a  reasonable  vertical  structure.  Further  indica¬ 
tion  of  their  similarity  is  contained  in  table  2.  The  ’’pillow"  referred 
to  is 

Pillow  (A  ~  (16) 

117  7 

and  the  RMSE  is  the  root-mean-square  error  calculated  in  the  usual 
fashion : 
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RMSE 


(17) 


rv>t  l*lT7  p  « 

P»IIonaJ 


\°\iy 


Nevertheless,  a  level-by- level  comparison  of  the  three  forecast 
methods  shows  some  differences  among  them.  At  1000  mb,  for  instance,  the 
two  low  centers  over  North  America  merged  and  deepened,  while  the  trough 
moved  nearly  20°  east.  All  methods  forecast  the  merging  and  consider¬ 
able  trough  movement,  with  method  I  giving  the  best  placement  of  the 
trough.  During  the  forecast  period  a  low  center  developed  in  the  area 
of  weak  gradient  southwest  of  the  Kamchatka  peninsula.  This  development 
was  forecast  by  all  methods,  although  with  too  much  cyclogenesis. 

Method  I  was  again  the  best.  In  no  case  is  there  any  evidence  of  spuri¬ 
ous  anticyclogenesis,  even  at  low  latitudes. 

The  forecasts  at  850  and  700  mb  show  movement  consistent  with  that 
forecast  at  1000  mb.  By  700  mb  the  tilt  of  the  baroclinic  systems  has 
become  apparent.  A  good  example  is  the  stationary  low  center  over  the 
Northwest  Territories,  which  deepened  by  26  meters.  All  forecasts  held 
it  stationary,  but  methods  I  and  III  forecast  two  and  seven  meters  deep¬ 
ening,  respectively,  while  method  II  forecast  a  2-meter  filling.  The 
weakening  of  the  low  center  east  of  Greenland  was  forecast  by  all  methods, 
but  best  by  method  II. 

At  500  mb,  since  the  chart  of  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  barotropic  forecast  has  been  included,  there  are  four  charts  to 
compare.  There  is  a  striking  similarity  between  the  barotropic  progno¬ 
sis  and  that  from  method  I.  The  method  II  and  III  prognoses,  themselves 
much  alike,  are  quite  different  from  the  other  two.  Method  I  and  the 
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USNFNWF  model  both  forecast  the  movement  and  reorientation  of  the  trough 
over  Eastern  United  States >  and  the  weakening  and  displacement  of  the 
low  center  near  Hudson  Bay*  Methods  II  and  III  show  the  trough  lagging 
far  behind  its  true  position,  and  retain  the  low  as  a  stationary  system* 
Method  III  even  deepened  it  slightly.  In  contrast,  methods  II  and  III 
produced  much  better  forecasts  than  the  others  in  the  western  Europe- 
Mediterranean  area. 

At  300  mb  all  three  methods  agree  in  general  on  the  movement  and 
development  of  the  major  systems.  This  development  is  in  some  cases  out 
of  phase  with  the  actual  development,  but  no  one  method  can  be  selected 
as  better  than  the  rest.  The  initial  and  verifying  charts  at  300  mb 
show  more  effect  from  the  process  of  making  the  data  vertically  consist¬ 
ent  than  was  evident  at  lower  levels. 

Referring  again  to  table  2,  there  seems  to  be  no  best  method  among 
the  three  tested.  Certainly  the  addition  of  the  advection  of  absolute 
vorticlty  by  the  divergent  wind  and  the  use  of  absolute  vorticity  in  the 
divergence  term  produced  little  change  and  no  over -ail  improvement  in 
this  particular  forecast*  The  effect  of  the  inclusion  of  the  vertical 
advection  of  absolute  vorticity  and  the  twisting  term  remains  to  be 
tested  * 
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Figure  1 


Table  1 


Level 

Contour  origin 

Contour  interval 

(mb) 

(meters) 

(meters) 

1000 

120 

60 

850 

1440 

60 

700 

3060 

60 

500 

5580 

60 

300 

9120 

120 

12 


Table  2 


Method 

Level 

(mb) 

Pillow 

(meters) 

msE 

(meters) 

I 

1000 

-  7.9 

43.6 

850 

-  1.8 

37.6 

700 

-  1.8 

41.8 

500 

-  1.2 

51.3 

300 

13.4 

80.3 

II 

1000 

-11.3 

43.3 

850 

-  2.7 

40.8 

700 

0 

42.7 

500 

-  .  6 

53.4 

300 

11.9 

77.3 

III 

1000 

-  9.1 

49.1 

850 

( 

-  3.4 

43.0 

700 

-  2.4 

43.6 

500 

-  1.5 

54.4 

300 

13.4 

80.0 

USNFNWF 

500 

-  4.9 

47.0 

Persistence 

500 

-  3.3 

57.4 
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APPENDIX  I 


The  following  charts  were  computed  during  the  course  of  this  study 
They  are  included  here  for  easy  reference. 


Charts  1-5 

Analyses ,  12Z  22  Nov  1962 

Charts  6-10 

24-hour  prognoses,  method  I 

Charts  11-15 

24-hour  prognoses,  method  II 

Charts  16-20 

24 -hour  prognoses 3  method  III 

Charts  21-25 

Analyses ,  122  23  Nov  1962 

Chart  26 

24-hour  prognosis,  USNFNWF 

15 


12Z-  22  NOU  C2 


CHART 


0?*8  * 

PROJECTION  POLAR  STEREOGRAPH  1C-  TRUE  AT  W  NORTH  LATITUDE 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTFRFV,  CALIFORNIA  ^ 


62  CHART  4 
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PROJECTION  POLAR  STEREOQRAPMIC-TMUE  AT  60'  >«ORTM  LATlTUOE 
SCALE  1  -60.000.000_ V 
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PROJECTION  POLAR  STEREOGRAPHIC- TRUE  AT  60"  NORTH  LATITUDE 
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PROJECTION  POLAR  STEREOGRAPHIC  -  TRUE  AT  60*  NORTH  LATITUDE 
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